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Abstract 

In this work, we develop and analyze a Hybrid High-Order (HHO) method for steady non¬ 
linear Leray-Lions problems. The proposed method has several assets, including the support for 
arbitrary approximation orders and general polytopal meshes. This is achieved by combining two 
key ingredients devised at the local level: a gradient reconstruction and a high-order stabilization 
term that generalizes the one originally introduced in the linear case. The convergence analysis is 
carried out using a compactness technique. Extending this technique to HHO methods has prompted 
us to develop a set of discrete functional analysis tools whose interest goes beyond the specific problem 
and method addressed in this work: (direct and) reverse Lebesgue and Sobolev embeddings for local 
polynomial spaces, L^-stability and kU^’^-approximation properties for L^-projectors on such spaces, 
and Sobolev embeddings for hybrid polynomial spaces. Numerical tests are presented to validate the 
theoretical results for the original method and variants thereof. 
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1 Introduction 


We are interested here in the numerical approximation of the steady Leray-Lions equation 

— div(a(-, M, Vw)) = / in O, (1.1a) 

u = 0 on do, 

where O c d ^ 1, is a polytopal bounded connected domain of boundary dO, while a : O x K x ^ 
is a (possibly nonlinear) function of its arguments, for which detailed assumptions are discussed in 


the following section. The homogeneous Dirichlet boundary condition (1.1b I is considered only for the 


sake of simplicity (the modifications required to handle more general boundary conditions are briefly 
addressed in the manuscript). This equation, which contains the p-Laplace equation, appears in the 
modelling of glacier motion [46| , of incompressible turbulent flows in porous media [35| a nd in airfoil 
design j^. Our goal is to design and analyze a discretization method for problem (1.1) inspired by 
the Hybrid High-Order (HHO) method introduced in 33 in the context of a linear diffusion model 
problem (see also for degenerate advection-diffusion-reaction models). The proposed method offers 
several assets:(i) the construction is dimension-independent; (ii) fairly general meshes including polytopal 
elements and nonmatching interfaces are supported; (iii) arbitrary polynomials orders can be considered 
(including the case k = 0); (iv) it is efficiently parallelisable (the local stencil only connects a mesh 
element with its faces), and it has reduced computational cost (when solving by a first-order algorithm, 
the element-based unknowns can be eliminated by static condensation). 
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^ jerome.droniou@monash.edu 
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Numerical methods allowing for arbitrary-order discretizations and general meshes have received 
increasing attention over the last few years. Supporting general polytopal meshes is required, e.g., in 
the modelling of underground flows, where degenerate elements and nonconforming interfaces account 
for complex geometric features resulting from compaction, erosion, and the onset of fractures or faults. 
Another relevant application of polyhedral meshes is adaptive mesh coarsening [^ |11| . The literature on 
arbitrary-order polytopal methods for linear diffusion problems is vast. In this context, methods that 
have similarities (and differences) with the HHO method include, e.g., the Hybridizable Discontinuous 
Galerkin method of (cf. also for a precise study of its relation with the HHO method), 

the Virtual Element Method of [T^[T^p^ , the High-Order Mimetic method of 51 
method of 54 55 , and the Multiscale Hybrid-Mixed method of [^. 


the Weak Galerkin 


The finite element approximation of nonlinear diffusion problems of Leray-Lions type on standard 
meshes has been studied in several papers; cf., e.g, 
is, 


10,46,52 


The literature on polytopal meshes 
however, much more scarce, and is mainly restricted to the lowest-order case. We cite here, in 


particular, the two-dimensional Discrete Duality Finite Volume schemes studied in (cf. also the 
precursor papers Bi). the Mixed Finite Volume scheme of (inspired by [^) valid in arbitrary 
space dimension, and the Mimetic Finite Difference method of for p G (1,2) and under more restrictive 
assumptions than ( |2.2[ ). High-order discontinuous Galerkin approximations have also been considered 
in [^. 


The starting point for the present work is the HHO method of 33 . In the lowest-order case, it has 


been shown in 33 Section 2.5] that this method belongs to the Hybrid Mixed Mimetic family 40 , which 
includes the mixed-hybrid Mimetic Finite Differences 20 , the Hybrid Finite Volume 43 and the Mixed 
Finite Volume [^. The HHO method can therefore be seen as a higher order version of these schemes. 
The (hybrid) degrees of freedom (DOFs) for the HHO method are fully discontinuous polynomials of 
degree fc > 0 at mesh elements and faces. The construction hinges on two key ingredients built element¬ 
wise: (i) a discrete gradient defined from element- and face-based DOFs; (ii) a high-order penalty term 
which vanishes whenever one of its arguments is a polynomial of degree ^ (fc -f 1) inside the element. 
These ingredients are combined to build a local contribution, which is then assembled element-wise. A 
key feature reducing the computational cost is that only face-based DOFs are globally coupled, whereas 
element-based DOFs can be locally eliminated by a standard static condensation procedure. 

The design of a HHO method for the nonlinear problem (1.1) entails several new ideas. A first 


difference with respect to the linear case is that a more natural choice is to seek the gradient reconstruction 
in the full space of vector-valued polynomials of degree ^ k (as opposed to the space spanned by gradients 
of scalar-valued polynomials of degree ^ {k + 1)). The main consequence of this choice is that, when 
applied to the interpolates of smooth functions, the discrete gradient operator commutes with the L^- 
projector, and therefore enjoys L^’-stability properties (see below). A second important point is the 
design of a high-order stabilization term with appropriate scaling. Here, we propose a generalization of 
the stabilization term of [33] which preserves the property of vanishing whenever one of its arguments 
is a polynomial of degree ^ (fc -f 1). As in the linear case, the construction hinges on the solution of 
small local linear problems inside each elements, and the possibility of statically condense element-based 
DOFs remains available. 


The convergence analysis is carried out using a compactness argument in the spirit of 53 . This 


technique, while not delivering an estimate of the convergence rate, has the crucial advantage of relying 
solely on the solution regularity inherent to the weak formulation. This point is particularly relevant 
for nonlinear problems, where additional regularity assumptions may turn out to be fictitious. The 
theoretical study of the convergence rate for smooth solutions is postponed to a future work. 

Adapting the compactness argument has prompted us to develop discrete functional analysis tools 
whose interest goes beyond the specific method and problem considered in this work. A first notable set 
of results are (direct and) reverse Lebesgue and Sobolev embeddings on local polynomial spaces (e.g., 
on mesh elements and faces, but curved geometries are also allowed). The term reverse refers to the fact 
that the largest exponent (semi-)norm is bounded above by the lowest exponent (semi-)norm. Direct 
Sobolev embedding for broken spaces on fairly general polytopal meshes are proved in 21 31 ; specific 
instances had already been established in [^[^|^[^[^. Reverse embeddings, on the other hand, are 


established in 18 Theorem 4.5.11], but under the assumption that all mesh elements are affine-equivalent 
to one (or a finite number of) given fixed reference elements. This limitation is due to the very generic 
local finite element spaces considered therein. Exploiting the fact that we deal with polynomial local 
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spaces, we can establish a more general version of reverse inequalities, that does not require to specify 
any particular geometry of the elements (only their non-degeneracy). Reverse Lebesgue embeddings are 
a crucial ingredient to prove the stability of the HHO method. 

A second set of results concerns the stability and approximation properties of the L^-projector on 
local polynomial spaces. More specihcally, we prove under very general geometric assumptions that the 
L^-projector is L^-stable for any index p e [1, -foo], and that it has optimal approximation properties in 
local polynomial spaces. Stability results for (global) projectors onto hnite element spaces can be found 
in iimii. However, these references mostly consider i/^-stability, and assume quite restrictive (and 
sometimes difficult to check) geometrical assumptions on the meshes. These limitations are a consequence 
of dealing with projectors on global hnite element spaces, that include some form of continuity property 
between the mesh elements. On discontinuous polynomial spaces such as the ones used in HHO methods, 
we can establish more general L^- and lT'’’P-stabihty and approximation properties of local L^-projectors. 


The approximation results extend to the lT®’^-setting the ones in 32 Section 1.4.4], based in turn on 


the ideas of 42 


Finally, a third set of discrete functional analysis tools are specihc to polynomial spaces with a hybrid 
structure, i.e., using as DOFs polynomials at elements and faces. In this case, building on the results 
for discontinuous Galerkin methods (inspired by the low-order discrete functional analysis results 


43 ), we introduce a suitable discrete hF^’^-like norm and prove a discrete counterpart of Sobolev 


of 
of 

embeddings and a compactness result for the discrete gradient reconstruction upon which the HHO 
method hinges. 

The material is organized as follows: in Section!^ we recall a set of standard assumptions to write 
a weak formulation for problem O; in Section we detail the discrete setting by specifying the 
assumptions on the mesh and recalling the basic results on local polynomial spaces; in Section we 
formulate the HHO method, state (without proof) the main stability and convergence results, and 
provide a few numerical examples; Section collects the discrete functional analysis tools on hybrid 
polynomial spaces, which are used in Section to prove the stability and convergence of the HHO 
method; in Section we briehy address the treatment of other boundary conditions and hint at the 
modihcations required in the analysis; a conclusion is given in Section and, hnally, in Appendix we 
provide the proofs of the discrete functional analysis results on local polynomial spaces. 


2 Continuous setting 

In this section we detail the assumptions on the function a and write a weak formulation for problem 
Let p e (1, -I-go) be given, and denote by p' := the dual exponent ofp, and by p* the Sobolev exponent 
of p such that 


I -1-00 if p > d. 


We assume that 


a : H X 


is a Caratheodory function. 


3ae LP'(H), e (0,-hoo), 3r < ^ : |a(a;, s, |)| ^ a(a:)-h +/3a|l|^ ^ 

for a.e. x e fl, for all (s, ^) e K x 

[sl{x, s, — a(a;, s, rj)] • [^ — J?] > 0 for a.e. x e fl, for all (s, rj) e K x x 

3Aa G (0, -l-oo) : a{x, s,$) ■ ^ ^ X e Q, for all (s, ^) e M x 

/eLP'(H). 


( 2 . 1 ) 

(2.2a) 

(2.2b) 

(2.2c) 

(2.2d) 

(2.2e) 


Here, Carathedory function means that a{x, •, •) is continuous on K x for a.e. a; e H, and a(-, s, is 
measurable on H for all (s, ^) e K x The Euclidean dot product and norm in are denoted by a; • y 
and I a; I, respectively. Classically 50 , the weak formulation for (1.1) is 


Find u G Wg^’^(H) such that, for ah v G Wq ’^(H), 
a{x,u{x),'Vu{x)) ■ Vv{x) dx = f{x)v{x)dx. 

Jn Jn 


(2.3) 
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The p-Laplace equation is probably the simplest type of Leray-Lions operator, and consists in setting 

a(a;, u, Vu) = I (2.4) 

In 14 , a simplified model of the stationary motion of glaciers is given by O) witii 

a{x,u,'Vu) = i^(|VM|)Vu, 


where F is the solution to the implicit equation F{s) ^ = {sF{s)) + Tq “ ; here, 


To > 0, and the unknown u in (1.1a) is the horizontal velocity of the ice. It is proved in 46 


a = 2^p G (0,1), 
that this 


choice of a satisfies (2.2). We refer the reader to 11^ for a discussion of models of turbulent flows using 


time-dependent versions of (1.1a) with a of the form 


a(a;, w, Vw) = |Vit —'^{Vu — h{u)) 


for some function ft, : K M'*. 

Existence of a solution to 
depend on s, the solution (whether weak or strong) is usually not unique, see e.g. 41, Remark 3.4]. 
Establishing a uniqueness result on (2.3) requires to strengthen the monotonicity assumption (2.2c). If 


2.3| is a consequence of the general results in 50 . Even if a does not 


a does not depend on s and is strictly monotone, in the sense that (2.2c) holds with a strict inequality 


whenever ^ =|= then the uniqueness of the solution to (2.3) is easy to see. Indeed, starting from two 


solutions u and u', subtracting the equations and taking v = u — u', we find 

[a(a;, VM(a;)) — a{x, VM'(a;))] • [VM(a;) — Vu'(a;)] dx = 0. 

Jo 

Since the integrand is non-negative, and strictly positive if 'Vu{x) Vu'(a;), this relation shows that 
Vm = Vu' a.e. on 11. We then deduce from the homogeneous boundary condition that u = u' a.e. on H. 
If a depends on s, the uniqueness of the solution is obtained by strengthening even more the monotonicity 


assumption (2.2c), and by assuming that a is Lipschitz continuous with respect to s, see 15 24 . 


3 Discrete setting 

This section presents the discrete setting: admissible mesh sequences, analysis tools on such meshes, 
DOEs, reduction maps, and reconstruction operators. 

3.1 Assumptions on the mesh 

Denote hy H cz a countable set of meshsizes having 0 as its unique accumulation point. Eollowing 
Chapter 4], we consider ft-refined mesh sequences {Th)he'H where, for all ft G "H, 71 is a finite collection 
of nonempty disjoint open polyhedral elements T such that fl = UreTh ^ ^ “ ^^^TeTh with 

ft-r standing for the diameter of the element T. A face F is defined as a hyperplanar closed connected 
subset of D with positive (d—l)-dimensional Hausdorff measure and such that (i) either there exist 
Ti , T 2 G 71 such that F cz dTi n 8X2 and F is called an interface or (ii) there exists T e Th such that 
F cz dT n dV, and F is called a boundary face. Interfaces are collected in the set boundary faces 
in and we let Fh '■= Fj^ F}^. The diameter of a face F e Fh is denoted hy hp. For all T G 71, 
Ft '■= {F e Fh \ F c dT} denotes the set of faces contained in dT (with dT denoting the boundary of 
T) and, for all F G Ft, riTF is the unit normal to F pointing out of T. Symmetrically, for all F e Fh, 
we let 7f := {T G 71 | T c: dT) the set of elements having T as a face. 

Our analysis hinges on the following assumption on the mesh sequence. 

Assumption 3.1 (Admissible mesh sequence). For all h e T-L, Fh admits a matching simplicial submesh 
T/i and there exists a real number g > 0 such that, for all ft G H: (i) for all simplex S B ^h of diameter 
hs and inradius rs, ghs ^ rs, and (ii) for all T G 71, and all S B Th such that S (z T, ghr < ft-s. 

The simplicial submesh in this assumption is just a theoretical tool, and it is not used in the actual 
construction of the discretization method. Given an admissible mesh sequence, for all ft G H, all T G 71, 
and all F B Ft, hp is uniformly comparable to hp in the sense that (cf. Lemma 1.42]): 

g^hp ^ hp ^ hp. (3-1) 


32 
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Moreover, 32, Lemma 1.41] shows that there exists an integer Ng depending on g such that 


Mhe'H ■. maxcardiJr) < Ng. 
TeTh 


(3.2) 


Finally, by 32, Lemma 1.40], there is an integer Ng depending on g such that 


'ihe'H : maxcard({S' e ] S' c T}) ^ Ng 
TeTh 


(3.3) 


3.2 Basic results on local polynomial spaces 

The building blocks for the HHO method are local polynomial spaces on elements and faces. Let an 
integer Z 5= 0 be fixed. Let 17 be a subset of (for some N > 1), Hjj the affine space spanned by U, 
du its dimension, and assume that U has a non-empty interior in Hjj. We denote by P*(?7) the space 
spanned by d( 7 -variate polynomials on Hu of total degree ^ Z. In the following sections, we will typically 
have N = d and the set U will represent a mesh element (and djj = cZ) or a mesh face (and djj = d — 1). 
We note, in passing, that a subset U with curved boundaries is also allowed except in Lemma |3.6[ which 
is why we use the different notation T instead of U in this lemma. 

A key element in the construction are L^-projectors onto local polynomial spaces on bounded subsets 
U c The L^-projector tt^ : L^{U) P*(Z7) is defined as follows: For any w e L^{U), is the 

unique element of P*(Z7) such that 

VueP*(Z7) : ( Tr[jw{x)v{x)dx = ( w{x)v{x)dx. (3.4) 

Ju Ju 

Note that the regularity w e L^(U) suffices to integrate w against polynomials on U (which are bounded 
functions). In what follows, we state some stability and approximation properties for the L^-projector. 
The proofs are postponed to Appendix |A.2[ 

Lemma 3.2 (LP-stability of L^-projectors on polynomial spaces). Let U be a measurable subset ofM.^, 
with inradius ru and diameter hjj, such that 


hu 


(3.5) 


Let fc e N and p g [1, -too]. Then, there exists C only depending on N, S, k and p such that 

ygeLP{U) : ^ C\\gU.iu)- 


(3.6) 

Remark 3.3 (Geometric regularity ( |3.5| ) for mesh elements and faces). Elements T e Th and faces 
F G Jpi of an admissible mesh sequence satisfy the geometric regularity assumption (3.51 with 5 = g^ and 
S = g respectively. 

In the case where IT®’P(Z7) is continuously embedded in C{U), the following result can be found 
1^ Theorem 4.4.4]. This restriction on the space IT®’P(Z7), which would prevent us from analyzing 


interesting cases for (1.11, is due to the very general setting chosen for analyzing the interpolation error. 
Because we focus here on local polynomial spaces and L^-projectors, we can improve this result and 
obtain optimal interpolation errors for any s,p. If U is an open set of s G N and p G [1, -too], we 
recall that j • |ws,p([/) is defined by 


Vu G IF*’P(Z7), 2 \Tv\\lhu), 

aeN-^, IckI^i =s 


where jaj^i = oi -I- ... -I- and d°‘ = ■ ■ ■ dff. 

Lemma 3.4 (IF®’P-approximation properties of L^-projectors on polynomial spaces). Let U be an open 
subset ofM.^ with diameter hjj, such that U is star-shaped with respect to a ball of radius phu for some 
p > 0. Let fc G N, s G {1,..., fc -I- 1} and p G [1, -too]. Then, there exists C only depending on N, p, k, s 
and p such that 

Vm G {0,..., s}, G W^’P(C/) : \v - ^ (3-7) 
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Figure 1: Degrees of freedom for k e {0,1, 2}. Shaded DOFs can be locally eliminated by static conden¬ 
sation. 


Remark 3.5. Using Section 7], the result still holds ifU is a finite union of domains that are star¬ 
shaped with respect to balls of radius comparable to hjj- This enables us to use Lemma M on elements 
of admissible mesh sequences, which are the union of a finite number of simplices; cf. (3.3). 

The next result estimates the trace of the error, and therefore requires more geometric assumptions 
on the domain (which, in the following sections, will be invariably a mesh element T). 

Lemma 3.6 (Approximation properties of traces of L^-projectors on polynomial spaces). Let T be a 
polyhedral subset ofMJ^ with diameter hx, such that T is the union of disjoint simplices S of diameter 
hs and inradius rs such that g^hx ^ ghs ^ rs for some g > 0. Let fc e N, s e {1,...,/c -t 1} and 
p e [1, -l-cx)]. Then, there exists C only depending on N, g, k, s and p such that 


Vtti s {0,..., s — 1}, Vu s (T) : hif\v — '^x^\w"'’‘i’(j^t) ^ Chfp '^\v\\Ys,ppxy 


(3.8) 


Here, VF'"’^(J^t) is the set of functions that belong to for any hyperplanar face F ofT, with 

corresponding broken norm. 


Finally, the triangle inequality applied to (3.7) (with m = s) and to (3.8) (with m = s—1) immediately 
gives the following extension of Lemma |3.2| 


Corollary 3.7 (bF^’^-stability of L^-projectors on polynomial spaces). The following holds: 

(i) Under the assumptions of Lemma \3.4\ we have, with C only depending on N, p, k, s and p, 

\lv£W’^’^{U) : \t^uv\w‘,v(u) ^ C\v\w=.p{u)\ 

(a) Under the assumptions of Lemma \3.()\ we have with C only depending on N, g, k, s andp, 
Vu e W^’P{T) : \TTxv\w=-^,p{rT) ^ klw-fCT) + 


4 The Hybrid High-Order method 

In this section we introduce the space of degrees of freedom, define the gradient and potential recon¬ 
structions at the heart of the HHO method, state the discrete problem along with the main stability and 
convergence results, and provide some numerical examples. 


4.1 Local degrees of freedom, interpolation and reconstructions 

Let a polynomial degree fc > 0 and an element T e Th he fixed. We define the local space of DOFs 

U^:=p'=(r)x| X P''(l^)), (4.1) 

cf. Figureand we use the underline notation Vx = (^t, (vf)f€J='t) ^ generic element y^ e U^. We 

define the local interpolation operator 1^ : W^-^{T) such that, for all v e W^'^{T), 

\^v := (tt^u, (7r^u)FeFT) ■ (4-2) 


6 











Remark 4.1 (Domain for the interpolation operator). The local interpolation operator is well-defined 
for functions v e since v is clearly in L^{T), the domain of irlf, and its trace on every face 

F G Ft is in L^{F), the domain ofiTp. In passing, in our convergence proofs we only need apply the 
interpolation operator to classically regular functions; cf, in particular, the proof of Theorem \4-()\ given 
in Section\^ 

Based on the local DOFs, we introduce reconstructions of the gradient and of the potential that will 
be instrumental in the formulation of the method. In what follows, (•, •)t and (•, ■)f denote the L^-inner 
products on T and F, respectively. The same notation is used in the vector case We define the 

local discrete gradient operator ^ such that, if := (vy, (vf)f€J='t) ^ then for all 

^gP'=(T)‘^, 


(G^y.j., c/))f — (Vvf, 0)t + ^ {\/f — ^T,(p-nTF)F 

= —(vf, V-^)^ + ^ {\/f,4>-'>T'Tf)i 
FeFt 


)F- 


(4.3a) 

(4.3b) 


Recalling the definition (4.2) of 1^, and using (4.3b I together with the definition ( |3.4[ ) of the L^-projector, 
one can prove that the following commuting property holds: For all v G VF^’^(T), 

Glfljfv = 4(Vr;), (4.4) 

where acts component-wise. As a result, by (3.7) and (3.8), Gp\Jf has optimal approximation prop¬ 
erties in P^(T)‘^. The local potential reconstruction operator : IJf ^ P^+^(T) is such that, for all 
Vp e U.T) gradient of is the orthogonal projection on VP^+^(T) of Gt)Lt^ and the average of 

p^+iy-p over T coincides with the average of vp, 

(Vp^'^^Vy — Gyy.j., Vic)t = 0 Vw g P*'''''^(T) and J {p^^\ip{x) —'jT{x))dx = Q. (4.5) 

For all V G H^{T), we have the following Euler equation: 

(V(p^+M^f-u),Vu;)t = 0 V'u;GP'=+i(r), (4.6) 

which shows that is nothing but the usual elliptic projector on P^+^(r). 


4.2 Global degrees of freedom, interpolation and reconstructions 

Local DOFs are collected in the following global space obtained by patching interface values: 

ufc j X P'=(r)) X ( X . 

\Ten / XFeFh / 

We use the notation y^ = {{\iT)TeTh^ {'^F)FeFh) ^ generic element y^ G and, for all T G Th, it 
is understood that Vp = (vf, (vf)f€J='t) denotes the restriction of y^ to U^. The global interpolation 
operator > ]J^ is defined such that, for all v G IF^’^(r2), 

((■^F'^lTeTfei ('^F'^)F€j='h)- (4-7) 

Interface DOFs are well-defined thanks to the regularity of functions in IF^d(r2). With V^{Th) usual 
broken polynomial space on Th, for all y^ G we denote by vu the unique function in P^(7)i) such that 

v/t|F = vf VT g Th- (4.8) 

Finally, we introduce the global discrete gradient operator Gfj : ^ ¥^{fThY ^•nd potential reconstruc¬ 
tion p^^ : ^ P^+^(7/i) such that, for all y^ G U^, 

(G^yJ,F = G^yr and {pl+\h)\T = p’F\t '^T^Th- (4.9) 
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4.3 Discrete problem and main results 


Define the following subspace of LJfi which strongly incorporates the homogeneous Dirichlet boundary 

(4.10) 


condition (1.1bI: 


Uto := y,eu^ 


\/F = 0 


We consider the following approximation of (2.3): 


Find U/j e U* Q such that, for any y^ G g, 4(u^,y;j) = f{x)vh{x)dx, (4.11a) 

Jo 


where 4 : U* x 


is assembled element-wise 


Ten 


(4.11b) 


from the local contributions At : IJt ^ Mt 


T e Th, defined such that 
At (Mt5 —t') ^ a(3?,U7^(ic), GrpUrp (x)) • Gt^Lt^^^^'^ ~b '^t(i1t'^ yr)? 

ST(u7’,yj.) := ^ f \T%{\iF — ^ 7r^(uF — ^’T''’^Uj.)(a;)7r^(vF — Py"''^yj.)(a;)ds(a;), 

VTci 7r„ F 


FeFt 


with : IJt ^ P^+^(T) denoting a second potential reconstruction such that, for all Vp ^ i=!.T> 


(4.11c) 


Pt'^^Vt ■= Vt + (PT^^yT ~ TpPt^^It)- 


(4.11d) 


Remark 4.2. This elaborate expression for the stabilization contribution sp aims at preserving the 
approximation qualities of the consistent contribution in Ap. As shown by ( |4.4[ ), Gp is exactly the 
gradient on (interpolations of) polynomials of degree ^ fc-f 1 inside the element. To preserve this exactness 
property in Ap, the stabilisation term sp must therefore vanish on (interpolations of) polynomials of 
degree ^ k + I inside the element. The choice in (4.11c) is one option that satisfies this property; 
other options include penalizing instead of Fp{\/F — PfA'^'Tp) o, combination of differences of the form 
Fpfvp ~ Pt^^)Lt) o-T-d Fp{\ip — p^^)Lt), weighted according the exponent p and their scaling properties 
with respect to the cell size. 

On the contrary, the more naive choice consisting in penalizing the difference (yp — vt) would only 
ensure that this stabilisation vanishes on polynomials of degree ^ k inside the element. This would 
prevent, e.g., from attaining the optimal convergence orders proved in l33fj for the linear case with p = 2. 


Remark 4.3 (Static condensation). Problem (4.11a) is a system of nonlinear algebraic equations, which 
can be solved using an iterative algorithm. When first order (Newton-like) algorithms are used, element- 
based DOFs can be locally eliminated at each iteration by a standard static condensation procedure. 

Remark 4.4 (Variants). Following 26) one could replace the space LJt of (4.1) with 


:= p«(r) X ) X 1P''(^) ( , 

for k ^ 0 and I e {k — l,k,k -\- 1}. For the sake of simplicity, we only consider here the case I = k — I 
when fc > 1. For k = 0 and I = k — 1, some technical modifications (not detailed here) are required 
owing to the absence of element-based DOFs. The local reconstruction operators Gp defined by (4.3) and 
Pp'^^ defined by (4.5) still map on P^(T)‘^ and P*"'"^(T), respectively (their domain changes, but we keep 
the same notation for the sake of simplicity). A close inspection shows that both key properties (4.4) 


and (4.6) remain valid for the proposed choices for 1. The second potential reconstruction operator 


defined by (4.lid), on the other hand, is replaced by 


l,k 


P^+i(T) such that, for all Vp e LJj. 


:= vp -f- (p^'^^Vp — TTpp^'^^v^p). The interest of the case I = k-\-l is that it holds, for allsup G Up, 


Pi 


,fc+l,fc+l 


)Lt = vt, o.nd the stabilization contribution takes the simpler form 


st{utNt) = S ^ f 


\p-2 


ttf^Uf - UT)(a:) 7r^(uF - UT)(a:)7r^(vF - yT){x)ds{x). 


FeFt 
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This simplification, however, comes at the price of having more element-based DOFs, which leads in 
turn to more onerous local problems for both the computation of the operator reconstructions and the 
elimination of element-based unknowns by static condensation. We also notice that the choice I = k + 1 
is close in spirit to the Hybridizable Discontinuous Galerkin methods introduced in for a linear 
diffusion problem. The choice I = k — 1, on the other hand, can be related to the High-Order Mimetic 
method introduced in [5^1/ in the context of linear elliptic equations. 


We next state our main results for problem (4.11). The proofs are postponed to Section]^ 


Theorem 4.5 (Existence of a discrete solution). Under Assumption (2.2), there exists at least one 


solution e Uf Q to (4.11). 


Theorem 4.6 (Convergence). We assume (2.2), and we let (Th)he'H be an admissible mesh seguence. 
For all h e H, we let e o be a solution to (4.11) on Th- Then up to a subsequence as h 0, 


recalling the definition (2.1) of the Sobolev index p*, 

• U/j —> It and —> u strongly in L'^{D) for all q < p*, 


• G'Ju 


h^h 


Vu weakly in 


where u e {Vl) solves the weak formulation (2.3) of the PDF (1.1). If we assume, moreover, that 


is strictly monotone, that is the inequality in (2.2c) is strict if ^ t], then 

Vu strongly in 




Remark 4.7 (Uniqueness). If a does not depend on s and is strictly monotone, then the solutions to 
both the continuous problem (2.3) and its discrete counterpart (4.11) are unique (see the discussion in 
Section^. In that case, the whole sequence of approximate solutions converges to the weak solution of 


Remark 4.8 (Other boundary conditions). The results stated in Theorems f.S-f.S are valid also when 


more general boundary conditions are considered (this is the case, e.g., in the numerical examples below). 
The modifications required to adapt the analysis to non-homogeneous Dirichlet and Neumann boundary 
conditions are briefly addressed in Section^ 


4.4 Numerical examples 

To close this section, we provide a few examples to numerically evaluate the convergence properties of 
the method (a theoretical study of the convergence rates is postponed to a future work). We consider the 


p-Laplace problem (2.4). When p = 2, we recover the usual (linear) Laplace operator, for which optimal 
convergence rates are proved in 33 . We consider the two-dimensional analytical solution originally 
proposed in Section 4], corresponding to u{x) = exp(xi 4- ttx 2 ) with suitable source term / inferred 
from ( l.la[). The domain is the unit square D = (0,1)^, and non-homogeneous Dirichlet boundary 


conditions inferred from the expression of u are enforced on its boundary; cf. (7.3) for the precise 


formulation of the method in this case. We compute the numerical solutions corresponding to polynomial 
degrees k = 0,..., 4. The meshes used are the triangular and Cartesian mesh families 1 and 2 from the 


FVCA 5 benchmark 47 , and the distorted (predominantly) hexagonal mesh family of 34, Section 4.2.3]; 
cf. Figure]^ 
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(a) Triangular mesh family (b) Cartesian mesh family (c) Hexagonal mesh family 


Figure 3: ||G^(u^ - [hu)\\Lp(ny vs. h, p = 3. 


-*-k = 0 -m-k=l -•-k = 2 -^k = 3 — fc = 4 





(a) Triangular mesh family (b) Cartesian mesh family (c) Hexagonal mesh family 

Figure 4: ||G^(u^ - lftM)||LP(n)d vs. h, p = 4. 


In Figures and 1^ we display the convergence of the error ||G*(u/j — i*M)||Lp(a)d for p = 3 and 
p = 4, respectively. In all the cases, we observe that increasing the polynomial degree k improves the 
convergence rate. The results obtained in 0il To] for lowest-order schemes suggest, however, that we 


should not expect optimal convergence properties in P^+^(7ii) except for the linear case p = 2. Instead, 
the order of convergence is expected to depend on both the regularity of the exact solution and the index 
p. Further numerical tests (not reported here for the sake of brevity) show that the convergence rate 
improves with k also when considering “degenerate” cases (i.e., solutions with a gradient that vanishes 
in part of the domain, in which case the diffusive properties of O) degenerate), although the gain is, in 
general, less relevant. Finally, for the sake of completeness, we report in Figure]^ the numerical results 
obtained for p = 4 with the method discussed in Remark 4.4 and corresponding to I = k + 1. In this 
case, taking the element-based DOFs in P^+^(r) does not seem to bring any significant advantage in 
terms of convergence (compare with Figure]^. 


5 Discrete functional analysis tools in hybrid polynomial spaces 

This section collects discrete functional analysis results on hybrid polynomial spaces that are used in the 
convergence analysis of Section 


5.1 Discrete IT^’^-norms 

We introduce the following discrete counterpart of the kF^’^-seminorm on ]jJ): 

hh\\i,p,h := f 2 > 

\Ten / 


(5.1) 
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(a) Triangular mesh family (b) Cartesian mesh family 



(c) Hexagonal mesh family 


Figures: ||G^(Uft - l^u)||Lp(n)<J vs. 
corresponding to ^ = A: + 1. 


h, p = A for the variant of the method discussed in Remark 


4.4 


and 


where the local seminorm on is defined by 




:= IIVvtI 


LP{T)‘ 




h]r^\\vF 


■ vt| 


LP{F) 


(5.2) 


It can be checked that the map defines a norm on IJ^q. We next show uniform equivalence 

between the local seminorm defined by (5.2) and two local kF^’P-seminorms defined using the discrete 
gradient and potential reconstructions (cf. ( |4.3a l and (4.5), respectively) and the penalty contribution 
St (cf. ( 4.11c[ )). This essentially proves stability for the discrete problem (4.11a) in terms of the 
norm. The argument hinges on the following direct and reverse Lebesgue embeddings, whose proof is 
postponed to Appendix |A.1| 

Lemma 5.1 (Direct and reverse Lebesgue embeddings). Let U be a measurable subset ofM.^ such that 
(3.5) holds. Let fc G N and g, m g [1,+oo]. Then, 


where A 
B sS MA. 


VruG ?'=([/) : \\wh.iu) ^ IC^Ilk|U™(c;), (5.3) 

B means that there is a real M > 0 only depending on N, k, S, q and m such that M~^A ^ 


We are now ready to prove the norm equivalence. 


Lemma 5.2 (Equivalence of discrete IT^’^-seminorms). Let {Th)he'H be an admissible mesh sequence 
and fc G N. Let T G 7h, P s [l,+oo), and denote by |•|s,p,T the local face seminorm such that, for all 
'It s U.T: recalling the definition (4.11c) of st, 


IVtIs.P.T := ST(yT,VT)^ = 2 h\,^\'K%{\JF 

XfsFt 


■ PT^y.T)i^)\^ds{x) 


(5.4) 


Then, 


tTU.P.T 


Vp‘, 


fe+1,, IIP 
'T 'tT\\LP(TY 


+ y7 


s.p.t) ^ (ll 


G^v, 


TXTllLP(T)'i 


+ V 7 


p y 

5,P,T J 


(5.5) 


where A ^ B means that M ^A ^ B ^ MA for some real number M > 0 that may depend on Ll, g, k 
andp, but does not otherwise depend on the mesh, T or\/rp. 


Remark 5.3 (Choice of the face seminorm). The proof of the norm equivalence does not make use of 
the specific structure of St, and could have been proved replacing Hs^p.T by any other local face seminorm 
composed by terms scaling on each face F e Ft as h^^^\\-\\Tp{F)- 
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Proof. We abridge A < B the inequality A ^ MB with real M only depending on g, k and p. 


Step 1: p = 2. It was proved in 33, Lemma 4] that 

I|VtIIi,2,T \\^PT^^y-T\\'L^{T)'i + \y-T\'i,2,Ti 


(5.6) 


which is exactly the first relation in (5.5) for p = 2. To prove the second, we notice that since, for 
all Vrp e LJy, VpIf^^Vj, is an orthogonal projection of in Lf(TY, we have \^P^^)Lt\\l‘^{tP ^ 

||Gyyj.||j;^ 2 (y)d. Relation (5.6) therefore shows that 

I|VtIIi, 2,T ~ l|G^TVTllL2(T)rf + llfTls,2,T- 


To prove the converse estimate, we make 4> = G^)tp into the definition (4.3a) of G^Vrp, and use the 
Cauchy-Schwarz inequality together with the discrete trace inequality 32 Lemma 1.46] to infer 

II^TlfTllL2(7’)d ^ II Vv7’j|^2(7n)d||G^yj.||^2(7n)d + ^ Ijv;’— Vr||i2(^) ||G^yj.||^2(7n)d 

F^J-t 

^ l|y'r||i,2,T||Gyy.p||2,2(r)d. 


This estimate shows that ||Gyy-p||/, 2 ( 7 ’)d < ||y 2 ’||i, 2 ,T and, combined with (5.6) to estimate |yy|s, 2 ,T 
||y 7 ’||i, 2 ,T, completes the proof of the case p = 2. 


< 


Step 2: p e [l,+oo). Relation (5.5) for a generic p can be deduced from the case p = 2 thanks to 
Lemma 5.1 (T and F clearly satisfy the geometric assumptions therein, cf. Remark 3.3). We only show 
how to do this to establish 


lll^Tlli,p,r ~ l|G'TVTllip('r)d + 


all the other estimates being obtained in a simil ar w ay. By admissibility of {Th)he'H^ we have hp\F\ ^ |T| 


for any F e Ft- Thus, for Vj, e ]J^, by Lemma |5.1[ 

l|vTlllp.T^|7’r'"l|VvT||^2(r)d+ 2 h]r^\F\^-^\vF-Ml^F) 


FeFt 


where, to pass to the second line, we used the inequality 

N 


^VT|li2(r)d + - Vr||i2(p.) 

F^Ft / 


N 


V6» > 0, Va* > 0 : Xl a* ^ 


a? 


(5.7) 


i=0 


which follows from writing aj = (a|)® < (Ziili ®i)® J- ^PP^y (5-5) with p = 2 and use again 

Lemma 5.1 and the inequality (5.7) to infer 


l|yTlll,p,T ~ ^ ( l|G^TyTllL2(7^)d + iT'F^W^pi'^F — Pt~''^'It)\\'l^(F) 


FeFt 


< 


\T\^ 


|T|^ p ||Gyyy||^p(j,)d + hp^\F\^ ’"W'^pi^F — PT'^^)tT)\\'LP{F) 


FeFt 


^ \T\ 


1 _£ 

o 


|T|’ 


+ Xj II'^f('^F PT~^^y.T)\\'LP(F) 


FeFt 


< ||G^y^||^p^.pjd + Xj ^F ^Kf(''f — £r''"^yT)llLP(F)- 
FeFt 


□ 
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5.2 Discrete Sobolev embeddings 


The first ingredient of our convergence analysis is the following discrete counterpart of Sobolev embed¬ 
dings, which will be used in Proposition |6.1| to obtain an a priori estimate of the discrete solution. 


Proposition 5.4 (Discrete Sobolev embeddings). Let {Th)hen be an admissible mesh sequence. Let 
\ ^ q ^ p* if 1 ^ p < d (with p* defined by (2.1)) and 1 ^ q < -l-oo if p ^ d. Then, there exists C only 
depending on LI, g, k, q and p such that 


Vy^ e 0 ■ l|v/i||L 9 (n) < 


(5.8) 


Remark 5.5 (Discrete Poincare). For q = p (this choice is always possible since p ^ p* for any space 
dimension d) this proposition states a discrete Poincare’s inequality. 

Proof. Here, A < B means that A ^ MB for some M only depending on LI, g, k, q and p. We recall the 
discrete Sobolev embeddings in P^(7/t) from 32 Theorem 5.3] (cf. also 21 31 ): 


ywe¥'^{Th) ■■ ||w||L<!(n) ^ ||w||dG.p, 

where the discrete W^’^-norm on P^(7)i) is defined by 


(5.9) 


llwjldG.p := Yi 
VTeTfe 




\\Lp{F) 


(5.10) 


Here, for all T e Th, wt '■= wy^, while [w]f •= wti — wt 2 is the jump of w through a face F e 
such that Tf = {Ti,T 2 } (the sign is irrelevant). If P e then Tf = {T} and we let [wJf = wt- For 
y^ e LJft 0 ^ ^ between Ti and T 2 , we have, using the triangle inequality, 

l|[Vh]F||LP(F) < I|vTi -Vf||fp(F) + II^Ta - Vf||lp(F)- 

Due to the strong boundary conditions, this estimate is also true if F is a boundary face and the term 
T 2 is removed. Hence, gathering by elements, 

E ^F”^IIK]F|lip(F) ^ E ^f"^ E llvT-VFllip(jp) = XI X ’T'F~^\\'^T-yF\\l^(^F)^hh\\ip,h- 

FeJ^h FeFh TeFf TeTh FeFt 


This shows that 


l|v/i||dG,p ^ l|y?illi,p,/i, 


which, plugged into (5.9), concludes the proof. 


(5.11) 

□ 


5.3 Compactness 

The second ingredient for our convergence analysis is the following compactness result for sequences 
bounded in the ll-jli^p^ft.-norm. 

Proposition 5.6 (Discrete compactness). Let {Fh)he'H be an admissible mesh sequence, and let, for all 
6 "H, y^ e LJ^ 0 be such that (||y/i||i.p,/t)teH is bounded. Then, there exists v e VFQ’^(n) such that, up to 
a subsequence as h^ 0, recalling the definition (2.1) of the Sobolev index p*, 


• '^h 


• C^y^ 


V and p\(^^)Lf, v strongly in L'^(Ll) for all q < p*, 
Vu weakly in LP{Lt)‘^. 


Remark 5.7. If p* < -l-oo, the discrete Sobolev embeddings (5.9) and Corollary 5.10 show that both V;, 
and p'F^y.h bounded in L^* (LI), and their convergence stated in Proposition 5.6 extends to L^* (Ll)- 
weak. 


The proof of Proposition 


5.6 


requires an auxiliary result allowing us to compare, for all y^ e LJfu the 


broken polynomial function (4.8) on Th defined by element DOFs and the potential reconstruction (4.9). 


Instrumental to obtaining this comparison result is the following Poincare-Wirtinger-Sobolev inequality 
on broken polynomial spaces, whose interest goes beyond the specific application considered here. 
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Lemma 5.8 (Poincare-Wirtinger-Sobolev inequality for broken polynomial functions with local zero 
average). Let {Th)hen be an admissible mesh sequence, and let p ^ q ^ p* with p* defined by (2.11. If 
w e P^(7)i) satisfies w{x)dx = 0 for all T e Th, then there exists C only depending on Ll, g, k, q and 
p such that (with V/i denoting the usual broken gradient), 


MlLiia) 






Lv(Q.Y 


(5.12) 


Remark 5.9. If p ^ d, the exponent 1 + ^ ~ ^ in h is positive if q < p* and equal to 0 if q = p*. 


have, for all T e Tk, tt^w = 0 and therefore, by (3.7) with fc = 0, s = 1 and m = 0, using Lemma 


Proof. In this proof, A < B means that A ^ AIB for some M only depending on O, g, k, q and p. We 

ir 

with m = p, and recalling that |T| < we write 


5.1 


1 _|_ “ _ “ 


\\w\\l'i{t) = \\w -'!Tfw\\Lq(^T) ^ hT\\'^w\\Lq(^T)'i ^ hT\T\i p j]V'u;||ip(T)d </ij, ® IIVr(;||ip(T)d- (5.13) 


If q is finite, we take the the power q of this inequality, sum over T e Th, and use || 
WVhwf^p^^y (we have q^ p) to infer 




w 


Li(n) 




Ten 


LP{T)‘ 




E l|V 

Ten 


W\ 


Lp{TY 






Vhwr-/^^y\\Vhw\%^^y = h^^'^- 




Vhw\\l,^^y. 


Taking the power 1/q of this inequality concludes the proof. If g = +oo, we apply (5.13) to T e Th such 
that ||r(;||icx)(T) = ||w||L«>(n) to obtain ||w||loo(q) < || Vw;||ip( 7 ^)d ^ || Vic|| 2 ,P(n)<J. □ 

Corollary 5.10 (Comparison between v^, and p^'^^Vff). Let {Th)he'H be an admissible mesh sequence, 
and let p q p*. Then, there exists C only depending on SI, g, k, q and p such that 


Vv/i e U,, : ||v?i— )Lh\\Li{n) ^ Ch “ p\\)Lh\\i,p,h- 


(5.14) 


Proof. Here, A < B means A ^ MB for M only depending on SI, g, k, q and p. By the second equation 
in (4.5), the average of Vh — Ph^^)Lh element of Th is zero. Hence, (5.12) gives 


l|v?* <h^+t ""^Whi^h- pt^^h^hpinY- 
Recalling the definitions (4.8) of V/, and of the ll-lli^p^ft.-norm, we have 


Ten 


(5.15) 


(5.16) 


Moreover, using the definition (4.9) of p^''’^y;, followed by the norm equivalence (5.5), and again the 
dehnition (5.1) of the H-Hi^p^/i-norm, it is inferred that 


2 < 2 llvrlllp.r = Kll?.p.,. 

Ten Ten 


(5.17) 


We conclude by using the triangle inequality in the right-hand side of (5.15) and plugging (5.16) and 
(5.17) into the resulting equation. □ 

We are now ready to prove the compactness result stated at the beginning of this section. 


Proof of Proposition \5.6[ By (5.11), (||v;i||dG,p)teW is bounded. The discrete Rellich-Kondrachov theo¬ 
rem Theor em 5.6] ensures that, up to a subsequence, Vh converges in L'?(SI) to some v. Since q < p*, 
Corollary |5 .1 0| shows that p^^y_h also converges in this space to the same v. 
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It remains to establish that v G VbQ’^(n) and that weakly converges to Vv. To this end, 

we first notice that G^y^ is bounded in U‘{Vt)'^ thanks to the norm equivalence (5.5). Hence, up to a 
subsequence, it weakly converges in to some Q. We take (/> G and observe that 


r Gf,yhi^yci>ix)dx = 
Jn 


= Y, (G^tYt - VVT, Cf> - 4c/))t + 2 (G^y^ - Vvt, 
Ten Ten 


+ 


2 (Vvt,</>)i 


Ten 


(cf. (E^) 


= Ti + ^ 2 {\/F —'^T,T'!y(f)-nTF)F + 2 (Vvr,(/))T 

Ten FeJ^T TeTh 

= h+ Y Xj (yp-'JT^i.^y^)-(j))-nTF)F - X (vT,div^)T (cf. ( |5.18[ )) 
Ten FePr Ten 


= Ti + T2 - 


v;i(a:) div <^(a:)(ia;. 

Jn 


In the penultimate line, we used a element-wise integration by parts, and the relation 

X X i^F,(i>-'>TTF)F = 0 , 

TeTh. FeJ^t 


(5.18) 


which follows from the homogeneous Dirichlet boundary condition incorporated in q (cf. (4.10)) and 
from riTiF + 'nT 2 F = 0 whenever F G is an interface between the two elements Ti and T 2 . If we prove 
that, as /i —» 0, Ti + T 2 ^ 0, then we can pass to the limit and we obtain 


Q{x) ■ 4>{x)dx = — v{x) dbf 4>{x)dx. 

Jn Jn 


(5.19) 


Taking (p compactly supported in H shows that G = Vn, and hence that v G IT^’P(H) and that G^y^ 
'Vv weakly in Taking then any cp G G®(K‘^)‘^ in ( |5.19[ ) and using an integration by parts shows 

that the trace of v on dVt vanishes, which establishes that v G Wq'^{GI). 

It therefore only remains to prove that Ti -f T 2 ^ 0. In what follows, A < B means that A ^ MB 
for some M not depending on h, cp or y^^. By Lemma 3.4 (with m = 0, s = 1 and p' instead of p) we 
have \(p — F^cplj^p!{ry ~ (Ty thus 


|Ti| </i ^ IG^y^-VvTlI 

\Ten 


P 

LP{Ty 


1/p 

I W'PWw'^-p'(ny ~ ^ {\\Gfi)Lh\\Lp{a)<i + ll<Mlwcp'(n)<*- 


Since ||y?i||i,p,?i is bounded, the norm equivalence (5.5) together with the definition (5.1) of the 
norm show that both ||Gj)y;j||ip(f 2 )d and || V?iV/j||j;,p(Q)d remain bounded. Hence, Ti 0 as /i ^ 0. The 
convergence analysis of T 2 is performed in a similar way. Using Lemma 3.6 (with p/ instead of p) we 

have \\<p — Lp'(F) ~ ^T’ll‘Mlivi'P'(T)-i thus, since /it ^ h-F whenever F G Ft, 


|‘I 2 | ^ X X ^F^F -'^t\\lp(F)\ 


(f>\\ 


wi.p'(r)‘ 


< 


TeTh FeTt 

~ M)^h\\l,P,h\\(p\\w'^’P' 

The convergence of T 2 to 0 follows. 


X X ^fI^f I ll 0 llwi,p'(a)<i 

\TeTh FeFt / 


□ 


5.4 Strong convergence of the interpolants 


The proof of Theorem 4.6 relies on a weak-strong convergence argument. The last ingredient of the 
convergence analysis is thus the strong convergence of both the discrete gradient and the stabilization 
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contribution when their argument is the interpolate of a smooth function. We state here this result in a 
framework covering more general cases than needed in the proof of Theorem 4.6 (where the argument of 
the interpolant is in C^(fi)). For r e N and g e [1, +oo], W^’'^(7h) denotes the broken space of functions 
(/? : n ^ M such that, for any T e 7h, iplr G This space is endowed with the norm 


M 


W^’i(Th) 




1/9 


\TeTh / 

miM 
Ih 


if g < +CX), 
if g = +CX). 


Proposition 5.11 (Strong convergence of interpolants). Let {Th)heH be an admissible mesh sequence, 
let p e [l,+oo], and let [f, be defined by (4.7). Then, there exists C not depending on h such that 


yg, e W^’\n) n W'^+^'P{Th) : - VpU.in) ^ 

As a consequence, 

yip G W^’P{n) : GI\Ip Vp strongly in LP{nY ash^Q. 

Moreover, 

yp G Wi’i(ff) n ■■ XI st{\^p,\Yp) ^0 ash 

TeTh 

Proof. We write A < B ior A ^ MB where M does not depend on h or p. 

Step 1: Proof of (5.20). By the commuting property (4.4) and the approximation property (3.7) 
applied to f = dip, s = k + 1 and m = 0, we have \\GY_tT ~ ^t\\lp{t)<^ ~ b,Y'^\\p\\w'‘+^’P(T) foi 
T e Th- Raising this inequality to the power p and summing over T e Tk (if p is finite, otherwise taking 
the maximum over T eTh) gives (5.20). 


(5.20) 

(5.21) 

(5.22) 


Step 2: Proof of (5.21). We reason by density. We take {pe)e>o c that converges to p in 

W^'P{Vl) as e —> 0 and we write, inserting ±{G^YiV<^ and using the triangle inequality. 


— V(p||ip(n)d sg \Gl\X{,p — Pe)||LP(n)d + IGlfifpe — ^Pe\\Lp(ny + l|V(P£ — p)llLp(n)'i 


-ik ik 


kik, 


< 


ll^(p ~ T<i)\Lp{aY + — Vipe||ip(n)d, 


where we have used the commuting property (4.4) followed by the Li’-stability of the L^-projector stated 
in Lemma 3.2 to pass to the second line. By (5.20), the second term in this right-hand side tends to 0 
as h —> 0. Taking (in that order) the supremum limit as /i —> 0 and then the supremum limit as e —> 0 
concludes the proof that G^Yp —> Vp in LP{LiY. 


Step 3: Proof of (5.22). It is proved in Eq. (46)] that 
hfYYFiO-TT)F - p^+%t)\\lhf) ^ 


Using Lemma 5.1 the admissibility of the mesh (which gives hp\F\ as |T| if Eg Et), and the regularity 


assumption on p, we infer 

\\Fp{{\xp) F — 


^\YFiilTT)F — PT^^[TT)\fLP(F) ^ b,p 


{hFYYF{i[TT)F - P^^YtT^LYF))' 


< {hF\F\y-^^hY^^YMi 


< \T\‘--ih. 


i-f f.(fc+i)p 
T 


Hk + 2(F) 

1^1 ^ l|p|l^fe + 2,oc.(7.) 


< 


|r|/r('=+i)p||pK 


Wfc+2,a>(r^)' 


Summing this inequality over F G Ft and T eTh, and recalling the uniform bound (3.2) over card(ET), 
we get 


X ST{\ip,YrP) 


Ten 


W'= + 2.aJ(rd)’ 


and the proof is complete. 


□ 
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6 Convergence analysis 


The following proposition contains an a priori estimate, uniform in h, on the solution to the discrete 


problem (4.111. 


3.1 


Proposition 6.1 (A priori estimates). Under Assumption 
exists C only depending on Vi, Aa, Q, k and p such that 




*/ llh e U^.o solves (|4.11[), then there 


( 6 . 1 ) 


Proof. We write A < B for A ^ MB with M having the same dependencies as C in the proposition. 


Plugging into (4.11a) and using the coercivity (2.2d) of a leads to 


TeTfi FeTt 


Ten 


Recalling the norm equivalence (5.5), and using the discrete Sobolev embeddings (5.8) with q = p to 


estimate the second factor in the right-hand side, this gives 

l|y/illi,p,/t ~ ll/llFp'(n)ll^'»IU^’(f2) ~ ll/llFp'(n)llilhlli.pA> 

which concludes the proof since, by assumption, p > 1. 


□ 


We can now prove that the discrete problem (4.11) has at least one solution. 


Proof of Theorem\4.5\ We use 29 Theorem 3.3] (see also 50): If (A, (•, •)f, H-Hf) is an Euclidean space. 


and ^ : E ^ E is continuous and satisfies ^ +oo as ||a;||F ^ +oo, then $ is onto. We take 

||1E|| p? 


E = LJ^ g, endowed with an arbitrary inner product, and define $ : LJ^ g —> LJ^ g by 

^)Lh,^h 6 U^_o, <d>(y^), w^>F = A(y^, w^). 


the norm equivalence (5.5) show that 


Assumptions (2.2a) and (2.2b) show that <i> is continuous, and the coercivity (2.2d) of a together with 


<$(yJ,y,>F > > CrJIyJI^, 

where Cji > 0 may depend on Th but does not depend on y^j (we use the equivalence of all norms on 
the finite-dimensional space LJ^ g). Hence, is onto. Let now e g be such that 

<y^, f{x)wh{x)dx Vw^ e g. 


and take e g such that <i>(u^) = y^. By definition of $ and y^, is a solution to the discrete 


problem (4.11). 


□ 


Let us now turn to the proof of convergence. To improve the legibility of certain formulas, we often 
drop the variable x inside integrals. 


Proof of Theorem\4.(H Step 1: Existence of a limit. By Propositions 


6.1 


and 


5.6 


there exists u e Wg’^(H) 

such that up to a subsequence as h ^ 0, Uh ^ u andp^''"^u,j ^ m in L'^(H) for all q < p*, and —> Vu 
weakly in Let us prove that u solves To this end, we adapt Minty’s technique to 

the discrete setting, as previously done in 36 41] . 


Step 2: Identification of the limit. The growth assumption (2.2b) on a ensures that a(-, u/j, G^u^) is 


bounded in (H)“, and converges therefore (upon extracting another subsequence) to some x weakly 
in this space. Let ip e G®(f2). Plugging y,j = \\ip into (4.11) gives 


r a(a;,u;„G^u,)-G^l^p= f st(u^,I^¥^) 

Jd 


( 6 . 2 ) 


Ten 
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with TT^ denoting the L^-projector on the broken polynomial space P^(7ft). Using Holder’s inequality 
followed by the norm equivalence (5.5) to bound the first factor, we infer 


Ten 


< 


ST(u7’,Ujn) j ( Xj ^t{[t^AtV^) 


\Ten 


\Ten 




y ^ stAtv’At^) 

\Ten j 


Recalling the a priori bound (6.1) on the exact solution and the strong convergence property (5.22), we 
see that this quantity tend s to 0 as /i —> 0. Additionally, by the approximation properties of the l?- 
projector stated in Lemma 3.4 together with the strong convergence property (5.21), we have 7r^(/3 ^ 
in LP(H) and ^ V(/3 in We can therefore pass to the limit ft, —> 0 in (6.2), and we find 


Jn J( 


(6.3) 


By density of in Wq’^(H), this relation still holds li ip £ WQ’^(r2). 

Let us now take A e L^iyiY write, using the monotonicity (2.2c) of a. 


I 


[a(a;, u/^, G^u^) - a(a;, Uh, A)] • [G^u^ - A] > 0. 


Use (4.11) and S 7 ’(up,Uj.) > 0 to write 


( a(a;, u?i, G^u,,) • G^u^ — | fuh 
Jn Jn 


Ten 


Develop (6.4) and plug this relation: 


A] > 0. 


(6.4) 

hr) ^ 

^ fUh. 

(6.5) 

J 



[gW ■ 

-A]. 

(6.6) 


Since u;, ^ m in L'^{Vl) for all q < p*, the Caratheodory and growth properties (2.2a) and (2.2b) of a 
show that a(a;,u/j. A) ^ a{x,u,A) strongly in We can therefore pass to the limit in (6.6): 


The conclusion then follows classically 


X • A 5= 


a(x, u, A) • [Vi 

Jn 


A]. 


(6.7) 


53 : Take v e IUq’^(U), apply this relation to A = Vm ± <Vu 
for some t > 0, use (|6.3|) with p = u±tv, divide by t, and let t ^ 0 using the Caratheodory and growth 


50 


properties of a. This leads to 


/u = a{x,u,Vu) 
Jn Jn 


Vu, 


and the proof that u solves (2.3) is complete. 

Step 3: Convergence of the gradient. It remains to show that if a is strictly monotone, then G^q^ 
Vw strongly in Let 


Fh = [a(x, u^, G^u^) - a(a;, u,,, Vu)] • [G^u^ - Vu] 3= 0 


( 6 . 8 ) 


Developing this expression and using (|6.5|), we can pass to the limit and use (6.3) to see that 

^u = 0. 


lim sup 

h^O 


f F/. < f /u - f X-Vz 

Jn Jn Jn 


Hence, F^ ^ 0 in Up to a subsequence, it therefore converges almost everywhere. Using the 

coercivity and growth assumptions (2.2d) and (2.2b) of a. Young’s inequality gives 

Ff, > AalG^uJP - (a(x) + + /3a|G^uJP-i)|Vu| - (a(a;) + + /3a|Vzr|P-i)|V m| 


s* ^\GWY - m^) + /3a|u,r)|Vzr| - /3a|Vur - - V 

2 p \P K 


p-i 


(6.9) 
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Since, up to a subsequence, converges a.e., this relation shows that for a.e. x, the sequence (G^u^(a;))/ig^ 
remains bounded. Let us show that it can only have Vu(a;) as adherence value. If C is an adherence 
value of then, passing to the limit in (| 6 . 8 |) gives, since Fh ^ 0 and Uh ^ u a.e.. 


[a(a;, u{x)X) — sl{x, u{x),'Vu{x))] ■ — Vu(a;)] = 0. 

The strict monotonicity of a then shows that C = Vu(x). Hence, for a.e. x, the bounded sequence 
(G^u^(x));ig'H has only Vu(a;) as adherence value, and thus G^u^ ^ Vu a.e. on fl. 

Since {Fh)h€'H is 1-equi-integrable (it converges in L^(n)) and {\uh\^)he'H is p'-equi-integrable {p'r < 
p* and {uh)he'H therefore converges in LP’'(n)), (6.9) shows that hen is p-equi-integrable. Vitali’s 

theorem then gives the strong convergence of this sequence to Vm in □ 


7 Other boundary conditions 

We briefly discuss here how the HHO scheme is written for non-homogeneous Dirichlet and homogeneous 
Neumann boundary conditions and hint at the modifications required in the convergence proof. 


7.1 Non-homogeneous Dirichlet boundary conditions 


Non-homogeneous Dirichlet boundary conditions consist in replacing (1.1b) with 

u = g on dn 


(7.1) 


with g e p’^{dQ). Denoting by 7 : p’^{dn) the trace operator, the weak formulation 

bGComGS ’ 

Find u e IF^’P(D) such that 7 ( 11 ) = g and, for all v e Wg^’*’(r2), 


a.{x,u{x),Vu{x)) ■'Vv{x)dx = f{x)v{x)dx. 

Jn Jn 


(7.2) 


As in Remark 4.1 we notice that -Kpg is well defined for any F e Fh. Hence, we can define the vector 


]Ag,h ^ Uft such that 


Ug.T = 0 VTeT^, Ug,F = 0 VFeJ-‘, = ^5 VF e 


We then set 


\lh,g \lh,0 + hg.fti 


and write the discrete problem corresponding to (7.2) as 


Find U;, e ^ such that, for any y^ e A(u^,y^) = fvh, 

Jn 


(7.3) 


with A defined by (4.11b)-(4.11c). The convergence analysis for non-homogeneous Dirichlet boundary 


conditions is performed as usual by utilizing a lifting of the boundary conditions. We take g e IF^’P(D) 


7.1 


and let = ^ 5 . Making y^ = e U;, in and using ||g^||i,p,?, < \\g\\w^’P{n) (see Proposition 


below) enables us to prove a priori estimates on ||u^ — 

Proposition |5.11| does not rely on the homogeneous boundary condi tions and therefore shows t hat 


Gpg —> V 5 in L^(ViY as h ^ Q. Since n^g 


to )Lh = 


g in LP(fl) (see Lemma 3.4), applying Proposition 
u 


5.6 


^ shows that, for some u e IT^’^(n) such that u — ge Wg^’^(fl) (i.e. 7 (u) = g), up to a 
subsequence u^, ^ m in L^’(H) and G^u^ ^ Vu in L^iVlY as /i ^ 0. The proof that u is a solution to 


(7.2) is then done in a similar way as for homogeneous boundary conditions. 


Proposition 7.1 (Discrete norm estimate for interpolate of IF^’^ functions). Let {Th)hen be an admis¬ 
sible mesh sequence, and let A: e N. Let v e and let e LJ^ be the interpolant defined by (4.7) 


and ( |4.2[ ). Then, ^ II^'IIivi.p(t) for all T eTh, and thus \\XAi,p,h ^ ll'f^llwi.p(a)• 
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LP(T) ^ 

Lemma 

13.21 with U = F and then Lemma 13.61 with m = 0 and s = 1 to write 

^ 

IIvf — vt||lp(f) = IKf'^ ~ ~ '^t'>^)\\lp(f) ^ 11^^ — '^t'^\\lp(f) ^ hj, ” ||r’||n/i,P(T)- 

Raising this to the power p and using hx ^ hp gives /if~^||vf — ^ ll'i'll^i,p(F)- global bound 

is then inferred raising the local bounds to the power p and summing over T e Fh- D 


Proof. Set v,j := \'fv and let T e Th- Since vp = f^v, Corollary 3.7 with s = 1 shows that [|Vvf|| 
Iklltvcp(T)- This takes care of the first term in p,T- To deal with the second term, we use 


7.2 Homogeneous Neumann boundary conditions 

We assume that 

f{x) dx = 0. 


I 

Jn 


Homogeneous Neumann boundary conditions for elliptic Leray-Lions problems consist in replacing (1.1b) 
with 

a(-, M, VM)-n = 0 on dn, (7.4) 


where n is the outer normal to dfl. The weak formulation of (l.la)-(7.4) is 

Find u e such that u{x) dx = 0 and, for all v e kF^’P(f2), 

Si{x,u{x),'Vu{x)) ■'Vv{x)dx = f{x)v{x)dx. 

Jn Jn 


(7.5) 


The HHO scheme for (7.5) reads 


Find G such that u?i(a;) da; = 0 and, for any G U^, H(u^, y^) = /v^ (7.6) 

Jn Jn 

with A still defined by (4.11b)-(4.11c). 

To carry out the convergence analysis from Section we need a few results. The first one is a 
discrete Poincare-Wirtinger-Sobolev inequality, which bounds to the -norm of discrete functions by 
their discrete norm. This immediately gives a priori estimates on the solution to the scheme (Proposition 


6.1). The second result is a discrete Rellich theorem for functions with zero average and bounded discrete 
norm (this is the equivalent of Proposition [5^. The proofs of both results are based on Lemma [5^ and 
on a decomposition of functions in LJ^ into low-order (piecewise-constant) vectors in and their higher 
order variation. 

Lemma 7.2 (Discrete Poincare-Wirtinger-Sobolev inequality for broken polynomial functions with zero 
global average). Let {Th)he'H be an admissible mesh sequence, and let q = p* if p A d, and q G [1,-l-cx)) 
if p = d. Then, there exists C only depending on Lt, g, k, q and p such that, for all Vf, e satisfying 

Vf,(x) dx = 0, we have 

l|v/t||L<!(n) < (7.7) 

Proof Here, A < B means that A sSi MB with M only depending on D, g, k and p. We define y° G ]J?1 
and vl G P^(Th) by: 

Vx = f^vt yTeTh, Vp = FpVp \/F e Fh, 

Vy = Vf — Fx'^x = Vf — Vf VT g Th- 


By Lemma 5.8 we have 


i/p 




< 


E iivvfI 

KTeTh 


P 

LP{T) 


(7.8) 


We recall the definition of the discrete W^’^-norm on LJ° from 


39 


hhWw^-p.u = ( Xj X 

XTeTh F€Ft 
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Assumption 


3.1 


(the genuine discrete W^^’^’-norm in 39 involves a different coefficient than |T| in this sum, but under 
this coefficient is « |r|). Since Yat^Th = ^Q'^h{x)dx = 0, 


39 


gives 


l|v° 11^9(0) SI IhhWw^-^Th- 


(7.9) 


By noticing that v/j = v); + the result follows from (7.81 and (7.9) provided that 


hhWw^-p.Th ~ l|V/illl,/i,p- 


(7.10) 


An easy generalisation of 37 Lemma 6.3] and 38 Lemma 6.6] (see 39 for details) shows that 

If If ^ hF 

| 7 r>T-<VT|^= j^^^^T{x)ds{x)-^^^^VT{x)dx < 1 ^ 

Using the triangular and Jensen’s inequalities, and the relations jTj < \F\hF and hp ^ hp, we infer 

hi 


J jVv7’(a;)j^da;. 


hP r 

- ^ J \^^T{x)\Pdx 


< 


1 r hP 

- Iv^(ai) - Mx)\^dsix) + ^11 


< 


^1 






Multiplying by ^ and summing over F e Ft and T e Fh gives (7.10). 


□ 


Proposition 7.3 (Compactness result for broken polynomial function with zero global average). Let 
('Th)he'H ^6 o,n admissible mesh sequence and let y^ e LJ^ be such that {\\'lh\\i,h,p)h.e'H 'Is bounded and, 
for all h € 'H, ^^Vh{x) dx = 0. Then, there exists v e W^’^{FV) such that ^^v{x)dx = 0 and, up to a 


subsequence as h —> 0, recalling the definition ( |2.1[ ) of the Sobolev index p*, 

• Vh —>■ V and p^'^^)/_h strongly in L'^(Cl) for all q < p*, 

• G^y^ ^ Vu weakly in LP{LI)‘^. 

Proof. We use the same decomposition as in the proof of Lemma 

we have I|v)jj|i9(f2) ^ G/i®j|y^j|i_ft_p where C does not depend on h and 6=1 + ^ ^ > 0. Hence, 

0 in L'^iFl) as h ^ 0. By (7.10), (l|y°l|wi.p,rfe)/te« remains bounded. Since 'ZaTbTh ^ 


7.2 


By Lemma 


5.8 


all /i e "H, the discrete compactness result for Neumann boundary conditions of shows that there 
exists a u G TU^’^(U) with zero average such that v° ^ v strongly in L'^(U) up to a subsequence. Hence, 
V/i ^ i; in L‘^(U) along the same subsequence. We then apply Corollary 5.10 which is independent of 
the boundary conditions, to deduce that p^^)Lh ^ in L'^{fl). 

To prove that G^y^^ 


Vu weakly in LP(HF, we notice that by Lemma 


5.2 


the functions 

remain bounded in LP{il)‘^ and ther efore converge weakly to some Q in this space. We prove that 
Q = Vu as in the proof of Proposition 5.6 using test functions (p e G®(U)‘’* instead of ^ G G”(M‘^)‘’*. □ 


8 Conclusion 


We extended the HHO method of 33 to fully non-linear Leray-Lions equations, which include the p- 


Laplace model. The lowest-order version of this method (corresponding to k = 0) belongs to the family 
of mixed-hybrid Mimetic Finite Differences, Hybrid Finite Volumes and Mixed Finite Volumes schemes. 
We proved the convergence of the HHO method without assuming unrealistic regularity properties on 
the solution, or restrictive assumptions on the non-linear operator. To establish this convergence, we 
developed discrete functional analysis results that include the analysis of L^- and IF®’^-stability and 
approximation properties of L^-projectors on broken polynomial spaces. We provided numerical results 
which demonstrate the good approximation properties of the method on a variety of meshes, and for 
various orders (low as well as high). 
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A Discrete functional analysis in local polynomial spaces 

This appendix collects discrete functional analysis results in local polynomial spaces that are of general 
interest for polynomial-based discretizations of linear and nonlinear problems. Most of these results have 
already been stated without proof in the paper, but we restate them for the sake of easy consultation. 


A.l Estimates in local polynomial spaces 

This section collects L^- and lT®’^’-estimates in local polynomial spaces including direct and reverse 
Sobolev and Lebesgue embeddings. 

Lemma 5.1 (Direct and reverse Lebesgue embeddings). Let U be a measurable subset ofM.^ such that 
(3.5) holds. Let k e N and q, me [1,+co]. Then, 


where A 
B sg MA. 


^we¥\U) : \\wU,iu) ^ \U\-^-^\\w\\l^^u), (5-3) 

B means that there is a real M > 0 only depending on N, k, 6, q and m such that M~^ A ^ 


Remark A.l (Reverse embeddings), //g ^ m then this result is a classical (direct) Lebesgue embedding 
due to Holder’s inequality. It holds for m < q solely because we consider polynomials (and we notice that 
the scaling explodes as hu 0). 

Remark A.2 (Sobolev reverse embeddings). Let U be a polyhedral set that admits a simplicial decom¬ 
position such that for any simplex S, if hs is the diameter of S and rg its inradius then hg ^ grg, and 
hu ^ ghg. The following inverse inequality holds with Cinv depending on g, k and p, but independent of 


h (cf. \3^ Lemma 1.44] for the case p = 2 and use use l3^ Lemma 1.50] or Lemma 5.1 to deduce the 
general case), 

\/neP^({7) : [|Vw||j;,p([/) ^ C'inv/lc/^||^'||Lp(c/)• (A.l) 

Using this inequality, we can easi ly d educe from Lemma \5.1\ the following reverse Sobolev embeddings: 
Under the assumptions of Lemma 5.1, if U is open and m ^ r, then for all w e P^({7) we have 

\w\w^’T{u) ^ h]f~'’’’^\U\p «|■^r|lyr■,<^(c/). 

Here < is up to a multiplicative constant only depending on k, 6, p, q and r. Note that the result obviously 
cannot hold if m < r and m ^ k (consider w polynomial of degree exactly m: the left-hand side does not 
vanish, while the right-hand side does). 


Proof of Lemma \5.1\ We obviously only have to prove < since m and q play symmetrical roles in (5.3). 
By (3.5), there is xu e U such that B{xij,Shij) cz U cz B(xu,hu). Let Uq = {U — xu)lhu. Using the 
change of variable x eU ^ {x — xij)/hij e U^, we see that, for i e [1, -l-oo], 




(A.2) 


where we used h(] « |17| (since hjj ss rjj) and we set WQ{y) = w{xij -f huv). Assume that there exists 
Co not depending on the geometry of Uq but solely on 6 such that 

VueP'"(t7o) : hWLiiUo) ^ ^o\\v\\l'^(Uo)- 


(A.3) 


Then combining this with (A.2), since wq e V^{Uo), 


I|w^IIl‘!(C/) ^ ’ l|li'o||L-!((7o) ~ \U\‘’‘\\'^o\\l^{Uo) ~ 1^1® ll^lk™((7)) 

and the lemma is proved. 

It remains to establish (A.3). To this end, we notice that, by choice of xjj, we have B{0,S) cz Co cz 
B{0, 1). Since ||•||L'J(B(o,l)) and |•||l™(b(o,i 5 )) are both norms on P^(?7o) (any polynomial that vanishes on 
a ball vanishes everywhere), and since P^([/o) is a finite-dimensional vector space, we have 

\/?;eP^(Co) II^IIl‘!(b(o,i)) ^ ll^^llz,’"(B(o,5))i (A.4) 


with constant in < depending on 6 but not on the geometry of Co. To prove (A.3), write 




□ 
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A.2 L^-stability and iy®’^-approximation properties of L^-projectors 

This section collects the proofs of L^- and IT^’^-stability and approximation estimates for L^-projectors 
on local polynomial spaces stated in Section [3^ 


Lemma 3.2 (LP-stability of L^-projectors on polynomial spaces). Let U be a measurable subset ofM.^, 
with inradius ru and diameter hjj, such that 


hu 


(3.5) 


Let fc e N and p e [1, +oo]. Then, there exists C only depending on N, 5, k and p such that 

yg G LP{U) : Wt^UWlhu) ^ C\\gU^uy (3-6) 

Proof. In this proof, A < B means that A ^ MB for some M only depending on N, S, k and p. 

Step 0: p = 2. This case is trivial since tt^ is an orthogonal projector in Lf{U) and therefore satisfies 
(3.6) with (7=1. 

Step 1: p > 2. We use Lemma 5.1 to write \\'n'^g\\Li‘(T) < \'r\p~^\\'^ug\\L^{T)- Since g e LP{T) c 

L^(T), we can use (3.6) for p = 2 and we deduce \\Tr^g\\Lp{T) ~ \T\p~^\\g\\L^(T)- We then conclude thanks 
to Holder’s inequality, valid since p > 2, 

\\^ug\\Lp(T) ^ \T\p ^\T\^ ” IIsIIlp(t) = llsllipcr)- 

Step 2: p < 2. We use a standard duality technique. Let g e LP{U) and w e L^ (U). Then by 
definition of tt^ and using (3.6) with p' > 2 instead of p, 

^^Tr^g{x)w{x)dx = ^^g{x)TT!fjw{x) dx ^ \\g\\Lp{u)huw\\Lp'(u) ^ \\g\\Lp(u)\\w\\Lp'(u)- 


Taking the supremum of this inequality over all w B L^ (U) such that ||'u;||= 1 shows that (3.6) 
holds. □ 

Lemma 3.4 (W®’P-approximation properties of L^-projectors on polynomial spaces). Let U be an open 
subset ofM.^ with diameter hjj, such that U is star-shaped with respect to a ball of radius phu for some 
p > 0. Let fc e N, s e {1,..., fc + 1} and p e [1, +cx)]. Then, there exists C only depending on N, p, k, s 
and p such that 


Vm e {0,..., s}, Vu e W^’P(C/) : \v - 7r^v\wpp.Piu) ^ Ch^^-^\v\ws.Piu)- 
Proof. Here, A < B means that A ^ MB with M only depending on N, p, k, s and p. 


(3.7) 



V = Q^v + 


(A.5) 


where Q^v is a polynomial of degree less than or equal to s— 1 and the remainder R^v satisfies 18, Lemma 
4.3.8] 

Vre{0,...,s} : {R^^ulw^^piu) ^ '^\v\w^.p{u)- (A-.6) 

Since is a polynomial of degree ^ s — 1 ^ fc, 7r^(Q^v) = Q^v and therefore, from (A.5), ir^v = 
Q^v + tt^{R‘^v). Subtracting this from (A.5), we infer v — tt^v = R^v — tt^{R‘^v). Hence, 


\v - TT^vlwp.piU) < |7?®u|w"'.p(( 7) + ky(i?®u)|w’".p(C/)- 


Iterating the inverse inequality (A.l) and using Lemma 3.2 we see that 


)(i?«u)| 


W'P’PiU) 


^ h. 


'riR%)\\ 


LP{U) ^ 


Estimate (A.6) applied to r = m and r = 0 shows that 

+ hf/^\R’'v\Lp{u) ^ h’'^™‘\v\w=’P{u)- 

The result follows from (A.7), ( |A.8 ) and (A.9). 


(A.7) 

(A.8) 

(A.9) 

□ 
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Lemma 3.6 (Approximation properties of traces of L^-projectors on polynomial spaces). Let T be a 
polyhedral subset of with diameter hx, such that T is the union of disjoint simplices S of diameter 
hs and inradius rs such that g^hx ^ ghs ^ rs for some g > 0. Let fc e N, s e {1,... ,A: + 1} and 
p e [1, +CX)]. Then, there exists C only depending on N, g, k, s and p such that 

Vme{0,...,s-1}, V?;GfL*’P(r) : h^\v - ^ Ch^^-"^\v\ws.P(x)- (3.8) 


Here, is the set of functions that belong to for any hyperplanar face F ofT, with 

corresponding broken norm. 


Proof. As expected A < is understood here up to a multiplicative constant that only depends on N, 
g, k, s and p. We first recall a classical continuous trace inequality: 


Vui e W^’^(T) : h^\\w\\xp{dT) < \\w\\lp{t) + hx\\^w\\xp{x)- 


(A.IO) 


For p = 2 this inequality can be deduced from 32 


Lemma 1.49] and many other references. The case 
of a general p is less easy to find in the literature, but actually very simple to prove. Since T is the 
union of disjoint simplices of inradius and diameter comparable to hx, it is sufficient to prove the result 
when T is one of these simplices S. For such a simplex, there exists an affine mapping A : T ^ Tq, 
where Tg = {a; e : Xi > 0, X;f=i < 1} i® reference simplex, such that the norms of the linear 
parts of A and A~^ are respectively of order and hx. Consider then wq e VF^’P(To) defined by 
wo(x) = w(A~^x). On To we have a trace inequality 


ll^«ollLP(eTo) ^ C!d,p{\\wo\\LP{Xo) + IIII L!>(To)) ■ 


(A.ll) 


By noticing that | Vu;o(®)| ^ hx\{Vw){A ^a;)| and using changes of variables x y = Ax, (A.ll ) gives 

dtrot . 

Estimate (3.8) is an immediate consequence of (A.IOI and of (3.7). For to ^ s —1, by applying (A.IO) 
to w = d°‘{v — TT^v) e IF^’P(T) for all a e of total length to we find 


h^\v - nxv\wrn,p{TT) ^ k - TTxi’lwp'-piT) + hx\v - 7 ry'u|^m+i,p(r). 


We then use (3.7) for to and to + 1 on the two terms in the right-hand side to conclude. 


□ 
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